gusucode.com > 现代通信系统——使用MATLAB（+全部程序） > 现代通信系统——使用MATLAB（+全部程序）/现代通信系统——使用MATLAB（+全部程序）/Matlab 程序/Chapter1/ip_01_02.m

    % MATLAB script for Illustrative Problem 2, Chapter 1.
echo on
fnct='lambda';
a=-4; 
b=4; 
n=24; 
tol=0.1;
xx=fseries(fnct,a,b,n,tol);
xx1=xx(n+1:-1:2);
xx1=[conj(xx1),xx];
absxx1=abs(xx1);
pause % Press any key to see a plot of the magnitude spectrum.
n1=[-n:n];
stem(n1,absxx1)
title('The Discrete Magnitude Spectrum')
phasexx1=angle(xx1);
pause % Press any key to see a plot of the phase.
stem(n1,phasexx1)
title('The Discrete Phase Spectrum')